This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue Howard Bayes
Time series 1999 - present 1977 - present
Bias in SWHS Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests. Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data. Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas. Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with. By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.
**Figure 1.**- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.


**Figure 2.**- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).


Results

**Figure X.**- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

**Figure 3.**- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 3.**- Total rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

**Figure 4.**- Black rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


And black rockfish releases…

**Figure 5.**- Black rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 6.**- Yellow rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 7.**- Yellow rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 8.**- DSR rockfish (excluding yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 9.**- DSR rockfish releases (including yelloweye) 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 11.**- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 12.**- Slope rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

**Figure 13.**- *Other* rockfish (non-pelagic and non-yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 14.**- *Other* rockfish (non-pelagic and non-yelloweye) releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



Total Biomass Removal Estimates

**Figure 15.**- Black rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.



**Figure 16.**- Yellow rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

**Figure 17.**- Pelagic rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 18.**- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 19.**- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


Model fit

Logbook residuals

**Figure 20.**- Residuals from logbook harvests.

Figure 20.- Residuals from logbook harvests.


SWHS residuals

**Figure 21.**- Residuals from SWHS harvests.

Figure 21.- Residuals from SWHS harvests.



**Figure 22.**- Residual of SWHS releases.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

**Figure 23.**- Mean percent of harvest by charter anglers.

Figure 23.- Mean percent of harvest by charter anglers.


When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

**Figure 24.**- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.


**Figure 25.**- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 26.**- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 27.**- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.


**Figure 28.**- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.


**Figure 29.**- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

**Figure 30.**- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.


Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

**Figure 31.**- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

**Figure 32.**- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

**Figure 33.**- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

**Figure 34.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

**Figure 35.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

**Figure 36.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.



P(slope|non-pelagic & non-yellowye) For release estimates

## NULL



Weight Fits

**Figure 37.**- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 38.**- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 39.**- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 40.**- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 41.**- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.4

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.05 while 0.95 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable

area_n

year

user

mean

median

lower

upper

Rhat

sd

skew

area

p_pelagic

13

1

1

0.322

0.309

0.105

0.614

1.042

0.135

0.096

NSEI

p_pelagic

13

1

2

0.140

0.125

0.036

0.338

1.053

0.079

0.189

NSEI

p_yellow

13

1

1

0.425

0.374

0.122

0.969

1.030

0.222

0.229

NSEI

p_yellow

13

1

2

0.430

0.379

0.136

0.968

1.039

0.219

0.233

NSEI

p_pelagic

13

10

1

0.336

0.318

0.118

0.627

1.051

0.136

0.132

NSEI

p_pelagic

13

10

2

0.146

0.130

0.041

0.345

1.054

0.080

0.199

NSEI

p_yellow

13

10

1

0.426

0.372

0.127

0.966

1.030

0.222

0.243

NSEI

p_yellow

13

10

2

0.428

0.377

0.129

0.971

1.035

0.223

0.229

NSEI

p_pelagic

13

11

1

0.342

0.330

0.120

0.636

1.045

0.134

0.090

NSEI

p_pelagic

13

11

2

0.147

0.130

0.040

0.348

1.043

0.080

0.213

NSEI

p_yellow

13

11

1

0.426

0.376

0.122

0.961

1.033

0.222

0.225

NSEI

p_yellow

13

11

2

0.427

0.378

0.134

0.974

1.031

0.222

0.221

NSEI

p_pelagic

13

12

1

0.345

0.329

0.123

0.636

1.034

0.136

0.117

NSEI

p_pelagic

13

12

2

0.148

0.131

0.040

0.339

1.048

0.080

0.213

NSEI

p_yellow

13

12

1

0.424

0.374

0.130

0.965

1.030

0.219

0.228

NSEI

p_yellow

13

12

2

0.426

0.379

0.121

0.969

1.033

0.221

0.213

NSEI

p_pelagic

13

13

1

0.353

0.345

0.124

0.636

1.033

0.134

0.060

NSEI

p_pelagic

13

13

2

0.149

0.133

0.044

0.354

1.046

0.081

0.198

NSEI

p_yellow

13

13

1

0.427

0.376

0.124

0.964

1.031

0.224

0.228

NSEI

p_yellow

13

13

2

0.431

0.383

0.130

0.970

1.034

0.221

0.217

NSEI

p_pelagic

13

14

1

0.330

0.317

0.128

0.608

1.033

0.124

0.105

NSEI

p_pelagic

13

14

2

0.148

0.132

0.046

0.339

1.040

0.078

0.206

NSEI

p_yellow

13

14

1

0.428

0.372

0.125

0.966

1.028

0.224

0.250

NSEI

p_yellow

13

14

2

0.429

0.379

0.126

0.968

1.030

0.221

0.226

NSEI

p_pelagic

13

15

1

0.305

0.288

0.113

0.581

1.034

0.123

0.138

NSEI

p_pelagic

13

15

2

0.142

0.129

0.044

0.320

1.035

0.073

0.179

NSEI

p_yellow

13

15

1

0.419

0.368

0.127

0.967

1.032

0.223

0.229

NSEI

p_yellow

13

15

2

0.425

0.374

0.125

0.969

1.033

0.221

0.230

NSEI

p_pelagic

13

16

1

0.344

0.331

0.132

0.627

1.027

0.128

0.101

NSEI

p_pelagic

13

16

2

0.154

0.141

0.047

0.337

1.035

0.077

0.169

NSEI

p_yellow

13

16

1

0.424

0.371

0.125

0.964

1.031

0.220

0.241

NSEI

p_yellow

13

16

2

0.428

0.374

0.128

0.969

1.034

0.222

0.243

NSEI

p_pelagic

13

17

1

0.353

0.342

0.146

0.619

1.026

0.122

0.090

NSEI

p_pelagic

13

17

2

0.161

0.147

0.052

0.349

1.026

0.078

0.180

NSEI

p_yellow

13

17

1

0.425

0.369

0.128

0.967

1.028

0.224

0.250

NSEI

p_yellow

13

17

2

0.430

0.376

0.130

0.970

1.031

0.222

0.243

NSEI

p_pelagic

13

18

1

0.362

0.353

0.147

0.617

1.033

0.124

0.073

NSEI

p_pelagic

13

18

2

0.165

0.149

0.053

0.364

1.024

0.080

0.199

NSEI

p_yellow

13

18

1

0.426

0.373

0.129

0.971

1.031

0.222

0.239

NSEI

p_yellow

13

18

2

0.430

0.380

0.124

0.972

1.025

0.222

0.225

NSEI

p_pelagic

13

19

1

0.360

0.349

0.157

0.622

1.022

0.123

0.089

NSEI

p_pelagic

13

19

2

0.166

0.152

0.054

0.356

1.026

0.078

0.178

NSEI

p_yellow

13

19

1

0.426

0.377

0.131

0.968

1.038

0.220

0.223

NSEI

p_yellow

13

19

2

0.428

0.374

0.124

0.968

1.032

0.221

0.245

NSEI

p_pelagic

13

2

1

0.324

0.305

0.108

0.638

1.053

0.139

0.137

NSEI

p_pelagic

13

2

2

0.138

0.120

0.036

0.340

1.061

0.080

0.226

NSEI

p_yellow

13

2

1

0.426

0.372

0.119

0.965

1.028

0.223

0.242

NSEI

p_yellow

13

2

2

0.427

0.379

0.128

0.970

1.033

0.221

0.218

NSEI

p_pelagic

13

20

1

0.372

0.361

0.158

0.625

1.025

0.122

0.090

NSEI

p_pelagic

13

20

2

0.173

0.162

0.058

0.365

1.018

0.079

0.139

NSEI

p_yellow

13

20

1

0.424

0.372

0.123

0.965

1.031

0.220

0.236

NSEI

p_yellow

13

20

2

0.428

0.373

0.128

0.968

1.034

0.218

0.252

NSEI

p_pelagic

13

21

1

0.389

0.380

0.180

0.645

1.019

0.122

0.074

NSEI

p_pelagic

13

21

2

0.183

0.168

0.064

0.377

1.018

0.082

0.183

NSEI

p_yellow

13

21

1

0.428

0.377

0.128

0.966

1.029

0.220

0.232

NSEI

p_yellow

13

21

2

0.427

0.375

0.124

0.968

1.028

0.222

0.234

NSEI

p_yellow

13

22

1

0.421

0.372

0.125

0.966

1.027

0.221

0.222

NSEI

p_yellow

13

22

2

0.431

0.386

0.131

0.969

1.037

0.219

0.206

NSEI

p_yellow

13

23

1

0.388

0.333

0.102

0.956

1.030

0.222

0.248

NSEI

p_yellow

13

23

2

0.480

0.439

0.154

0.977

1.026

0.218

0.188

NSEI

p_yellow

13

24

1

0.397

0.341

0.099

0.958

1.040

0.224

0.250

NSEI

p_yellow

13

24

2

0.473

0.435

0.159

0.975

1.028

0.214

0.177

NSEI

p_yellow

13

25

1

0.394

0.340

0.105

0.962

1.041

0.222

0.243

NSEI

p_yellow

13

25

2

0.471

0.430

0.159

0.972

1.028

0.216

0.190

NSEI

p_yellow

13

26

1

0.388

0.333

0.106

0.963

1.035

0.221

0.248

NSEI

p_yellow

13

26

2

0.476

0.438

0.153

0.972

1.023

0.216

0.176

NSEI

p_yellow

13

27

1

0.374

0.321

0.102

0.954

1.031

0.214

0.247

NSEI

p_yellow

13

27

2

0.468

0.431

0.152

0.972

1.025

0.212

0.175

NSEI

p_yellow

13

28

1

0.388

0.342

0.108

0.950

1.039

0.209

0.220

NSEI

p_yellow

13

28

2

0.468

0.434

0.160

0.969

1.026

0.208

0.163

NSEI

p_yellow

13

29

1

0.361

0.318

0.099

0.903

1.028

0.198

0.217

NSEI

p_yellow

13

29

2

0.459

0.433

0.157

0.942

1.026

0.198

0.131

NSEI

p_pelagic

13

3

1

0.328

0.309

0.110

0.639

1.050

0.139

0.137

NSEI

p_pelagic

13

3

2

0.140

0.122

0.036

0.340

1.039

0.079

0.227

NSEI

p_yellow

13

3

1

0.426

0.371

0.128

0.967

1.037

0.220

0.250

NSEI

p_yellow

13

3

2

0.434

0.380

0.136

0.966

1.029

0.219

0.247

NSEI

p_pelagic

13

4

1

0.328

0.313

0.111

0.633

1.056

0.136

0.110

NSEI

p_pelagic

13

4

2

0.140

0.122

0.036

0.339

1.044

0.080

0.225

NSEI

p_yellow

13

4

1

0.427

0.376

0.126

0.966

1.031

0.223

0.229

NSEI

p_yellow

13

4

2

0.431

0.382

0.125

0.968

1.033

0.220

0.223

NSEI

p_pelagic

13

5

1

0.331

0.314

0.116

0.624

1.051

0.136

0.125

NSEI

p_pelagic

13

5

2

0.143

0.125

0.037

0.344

1.052

0.081

0.221

NSEI

p_yellow

13

5

1

0.425

0.377

0.119

0.966

1.033

0.222

0.216

NSEI

p_yellow

13

5

2

0.430

0.381

0.127

0.969

1.028

0.224

0.219

NSEI

p_pelagic

13

6

1

0.330

0.313

0.109

0.637

1.038

0.139

0.123

NSEI

p_pelagic

13

6

2

0.144

0.127

0.037

0.343

1.048

0.082

0.207

NSEI

p_yellow

13

6

1

0.427

0.374

0.122

0.963

1.036

0.221

0.240

NSEI

p_yellow

13

6

2

0.426

0.377

0.123

0.966

1.030

0.224

0.219

NSEI

p_pelagic

13

7

1

0.335

0.320

0.119

0.629

1.037

0.137

0.109

NSEI

p_pelagic

13

7

2

0.143

0.125

0.038

0.343

1.046

0.080

0.225

NSEI

p_yellow

13

7

1

0.425

0.372

0.122

0.966

1.033

0.221

0.240

NSEI

p_yellow

13

7

2

0.429

0.376

0.126

0.971

1.030

0.224

0.237

NSEI

p_pelagic

13

8

1

0.331

0.317

0.115

0.621

1.044

0.136

0.103

NSEI

p_pelagic

13

8

2

0.146

0.127

0.038

0.354

1.047

0.082

0.232

NSEI

p_yellow

13

8

1

0.425

0.378

0.128

0.967

1.037

0.220

0.213

NSEI

p_yellow

13

8

2

0.429

0.378

0.135

0.969

1.036

0.219

0.233

NSEI

p_pelagic

13

9

1

0.336

0.321

0.115

0.632

1.036

0.137

0.109

NSEI

p_pelagic

13

9

2

0.146

0.127

0.039

0.345

1.044

0.081

0.234

NSEI

p_yellow

13

9

1

0.426

0.369

0.123

0.967

1.033

0.223

0.255

NSEI

p_yellow

13

9

2

0.433

0.381

0.129

0.970

1.027

0.223

0.233

NSEI

p_pelagic

16

1

1

0.375

0.369

0.134

0.665

1.010

0.140

0.043

SSEO

p_pelagic

16

1

2

0.349

0.338

0.125

0.633

1.013

0.134

0.082

SSEO

p_pelagic

16

10

1

0.378

0.369

0.150

0.663

1.013

0.136

0.066

SSEO

p_pelagic

16

10

2

0.355

0.344

0.129

0.633

1.016

0.134

0.082

SSEO

p_pelagic

16

11

1

0.380

0.370

0.147

0.662

1.011

0.135

0.074

SSEO

p_pelagic

16

11

2

0.355

0.344

0.129

0.646

1.018

0.135

0.082

SSEO

p_pelagic

16

12

1

0.379

0.368

0.152

0.649

1.016

0.132

0.083

SSEO

p_pelagic

16

13

2

0.364

0.353

0.135

0.638

1.014

0.132

0.083

SSEO

p_pelagic

16

15

1

0.375

0.367

0.155

0.636

1.011

0.128

0.063

SSEO

p_pelagic

16

2

2

0.352

0.343

0.119

0.634

1.020

0.134

0.067

SSEO

p_pelagic

16

3

1

0.372

0.364

0.137

0.651

1.014

0.135

0.059

SSEO

p_pelagic

16

3

2

0.352

0.343

0.123

0.625

1.013

0.134

0.067

SSEO

p_dsr

16

31

1

0.399

0.398

0.375

0.431

1.270

0.014

0.070

SSEO

p_slope

16

36

1

0.062

0.062

0.044

0.079

1.027

0.008

0.000

SSEO

p_pelagic

16

4

1

0.378

0.371

0.139

0.660

1.014

0.135

0.052

SSEO

p_pelagic

16

4

2

0.352

0.345

0.129

0.634

1.013

0.134

0.052

SSEO

p_pelagic

16

42

1

0.832

0.832

0.826

0.839

1.103

0.003

0.000

SSEO

p_yellow

16

42

1

0.405

0.405

0.383

0.427

1.040

0.011

0.000

SSEO

p_dsr

16

42

1

0.469

0.479

0.411

0.513

3.654

0.030

-0.332

SSEO

p_slope

16

42

1

0.146

0.133

0.106

0.205

4.921

0.031

0.415

SSEO

p_pelagic

16

5

1

0.371

0.364

0.143

0.646

1.012

0.131

0.053

SSEO

p_pelagic

16

5

2

0.354

0.344

0.121

0.631

1.022

0.136

0.074

SSEO

p_pelagic

16

6

1

0.372

0.365

0.137

0.653

1.013

0.136

0.052

SSEO

p_pelagic

16

6

2

0.351

0.342

0.126

0.641

1.012

0.135

0.067

SSEO

p_pelagic

16

7

1

0.377

0.370

0.141

0.662

1.011

0.134

0.052

SSEO

p_pelagic

16

7

2

0.355

0.344

0.123

0.641

1.017

0.134

0.082

SSEO

p_pelagic

16

8

1

0.373

0.364

0.139

0.661

1.010

0.135

0.067

SSEO

p_pelagic

16

8

2

0.354

0.345

0.127

0.638

1.023

0.136

0.066

SSEO

p_pelagic

16

9

1

0.377

0.370

0.144

0.662

1.014

0.135

0.052

SSEO

p_black

8

40

1

0.488

0.418

0.343

0.963

1.135

0.191

0.367

afognak

p_black

8

48

1

0.531

0.492

0.414

0.967

1.024

0.138

0.282

afognak

p_black

9

32

1

0.912

0.942

0.664

0.993

1.011

0.091

-0.330

eastside

Table Xb. True parameters with an Rhat in excess of 1.1

variable

area_n

year/species

mean

median

lower

upper

Rhat

sd

skew

year

area

beta2_pH

10

1

1.434

0.855

0.320

6.979

1.123

1.762

0.329

northeast

beta0_pelagic

13

-0.833

-0.812

-1.682

-0.001

1.110

0.464

-0.045

NSEI

beta3_pelagic

13

31.843

31.539

24.861

40.019

1.116

4.165

0.073

NSEI

beta3_pH

15

4

41.375

41.277

39.527

43.390

1.216

1.012

0.097

SSEI

beta3_pH

16

4

41.536

41.522

40.923

42.287

1.257

0.368

0.038

SSEO

beta0_dsr

16

0.165

0.164

0.133

0.196

1.150

0.016

0.062

SSEO

mu_beta2_pH

2

1

1.436

0.917

0.205

5.621

1.128

1.710

0.304

NG

beta2_pH

4

1

0.952

0.702

0.184

2.866

1.126

0.733

0.341

PWSO

beta1_pH

6

1

2.078

1.866

1.045

4.236

1.122

0.794

0.267

SOKO2SAP

beta2_pH

6

1

1.185

0.510

0.120

6.945

1.136

1.963

0.344

SOKO2SAP

beta2_pH

9

1

1.179

0.630

0.173

6.278

1.126

1.676

0.328

eastside

lambda_H

10

0.233

0.154

0.023

0.848

1.128

0.447

0.177

northeast

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.01

538.00

222.00

1,032.00

353.00

Prop. of total

0.71

0.29

0.75

0.25

Prop. median less than 1

0.27

0.22

0.34

0.42

Prop. median less than 10

0.39

0.55

0.55

0.73

Prop. median less than 50

0.51

0.79

0.70

0.82

Median of the median values

45.65

8.03

5.66

1.98

Mean of the median values

549.33

354.99

179.87

121.86

Mean skew of posteriors

0.25

0.29

0.27

0.25

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.1

7.00

4.00

70.00

36.00

Prop. of total

0.64

0.36

0.66

0.34

Prop. median less than 1

0.43

1.00

0.73

0.69

Prop. median less than 10

0.43

1.00

0.84

1.00

Prop. median less than 50

0.43

1.00

0.91

1.00

Median of the median values

130.65

0.10

0.03

0.09

Mean of the median values

609.22

0.13

33.47

1.05

Mean skew of posteriors

0.25

0.20

0.15

0.16

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable

Freq

beta0_dsr

1

beta0_pelagic

1

beta1_pH

1

beta2_pH

4

beta3_pelagic

1

beta3_pH

2

eps_bl

1

eps_dsr

1

eps_pH

8

eps_sl

2

H_ay

2

H_ayg

1

H_ayu

1

Hb_ay

2

Hb_ayg

1

Hb_ayu

1

Hd_ay

1

Hd_ayg

2

Ho_ay

4

Ho_ayu

4

Hp_ay

2

Hp_ayg

1

Hp_ayu

1

Hs_ay

2

Hs_ayg

2

Hs_ayu

2

Hy_ayu

1

lambda_H

1

mu_beta2_pH

1

p_black

1

p_dsr

2

p_slope

2

pH

14

R_ayg

1

R_ayu

3

Rb_ayg

2

Rd_ay

5

Rd_ayg

6

Rd_ayu

2

re_black

1

re_dsr

2

re_pH

11

re_slope

2

Ro_ay

16

Ro_ayg

11

Ro_ayu

43

Rp_ayg

1

Rs_ay

9

Rs_ayg

13

Rs_ayu

20

Ry_ay

12

Ry_ayg

8

Ry_ayu

26

Parameter estimates:

Summary Table of Parameter Estimates
Parameter mean sd Lower_CI Median Upper_CI
mu_bc_H[1] -0.019 0.069 -0.146 -0.022 0.124
mu_bc_H[2] 0.069 0.054 -0.031 0.068 0.180
mu_bc_H[3] -0.279 0.067 -0.407 -0.280 -0.142
mu_bc_H[4] -0.922 0.185 -1.286 -0.921 -0.572
mu_bc_H[5] 1.344 1.042 0.129 1.084 3.985
mu_bc_H[6] -2.375 0.319 -2.983 -2.381 -1.741
mu_bc_H[7] -0.455 0.101 -0.664 -0.451 -0.267
mu_bc_H[8] 0.253 0.339 -0.307 0.225 1.019
mu_bc_H[9] -0.245 0.123 -0.489 -0.247 0.004
mu_bc_H[10] 0.059 0.064 -0.060 0.058 0.186
mu_bc_H[11] -0.126 0.037 -0.201 -0.126 -0.054
mu_bc_H[12] -0.247 0.108 -0.479 -0.245 -0.053
mu_bc_H[13] -0.163 0.073 -0.302 -0.164 -0.017
mu_bc_H[14] -0.307 0.090 -0.493 -0.305 -0.133
mu_bc_H[15] -0.336 0.049 -0.432 -0.336 -0.242
mu_bc_H[16] -0.152 0.384 -0.794 -0.184 0.724
mu_bc_R[1] 1.044 0.120 0.801 1.046 1.271
mu_bc_R[2] 1.009 0.077 0.852 1.009 1.156
mu_bc_R[3] 0.236 0.093 0.048 0.238 0.416
mu_bc_R[4] 0.674 0.196 0.277 0.679 1.049
mu_bc_R[5] 1.681 0.355 0.982 1.685 2.360
mu_bc_R[6] -1.415 0.343 -2.062 -1.423 -0.718
mu_bc_R[7] 0.509 0.168 0.164 0.512 0.842
mu_bc_R[8] 0.652 0.143 0.374 0.650 0.943
mu_bc_R[9] 0.401 0.164 0.050 0.408 0.697
mu_bc_R[10] 0.902 0.102 0.697 0.904 1.099
mu_bc_R[11] 1.108 0.074 0.960 1.107 1.255
mu_bc_R[12] 0.881 0.161 0.567 0.880 1.200
mu_bc_R[13] 1.042 0.091 0.859 1.044 1.220
mu_bc_R[14] 0.867 0.121 0.624 0.867 1.095
mu_bc_R[15] 0.800 0.098 0.606 0.800 0.992
mu_bc_R[16] 1.359 0.115 1.133 1.359 1.582
tau_pH[1] 2.432 0.179 2.091 2.430 2.809
tau_pH[2] 1.635 0.157 1.336 1.635 1.950
tau_pH[3] 1.220 0.112 1.018 1.216 1.458
tau_pH[4] 1.109 0.149 0.847 1.102 1.430
tau_pH[5] 0.438 0.078 0.303 0.431 0.609
beta0_pH[1,1] 0.544 0.184 0.143 0.557 0.876
beta0_pH[2,1] 1.391 0.167 1.058 1.392 1.704
beta0_pH[3,1] 1.377 0.256 0.753 1.408 1.761
beta0_pH[4,1] 1.536 0.251 0.966 1.542 1.988
beta0_pH[5,1] -0.763 0.373 -1.614 -0.729 -0.130
beta0_pH[6,1] -0.515 0.479 -1.931 -0.388 0.065
beta0_pH[7,1] 0.759 0.011 0.739 0.759 0.780
beta0_pH[8,1] -0.510 0.300 -1.143 -0.495 0.006
beta0_pH[9,1] -0.305 0.269 -0.872 -0.287 0.181
beta0_pH[10,1] 0.280 0.196 -0.129 0.293 0.643
beta0_pH[11,1] 0.094 0.325 -0.563 0.116 0.672
beta0_pH[12,1] 0.609 0.191 0.213 0.610 0.970
beta0_pH[13,1] -0.655 0.151 -0.985 -0.645 -0.390
beta0_pH[14,1] 0.020 0.173 -0.326 0.022 0.347
beta0_pH[15,1] 0.047 0.198 -0.392 0.063 0.384
beta0_pH[16,1] -0.858 0.216 -1.311 -0.852 -0.453
beta0_pH[1,2] 2.565 0.084 2.408 2.563 2.728
beta0_pH[2,2] 3.079 0.022 3.037 3.078 3.122
beta0_pH[3,2] 3.626 0.041 3.549 3.625 3.709
beta0_pH[4,2] 2.936 0.019 2.898 2.936 2.975
beta0_pH[5,2] 4.518 1.462 2.761 4.223 7.960
beta0_pH[6,2] 3.326 0.156 3.028 3.325 3.650
beta0_pH[7,2] 2.034 0.051 1.936 2.033 2.133
beta0_pH[8,2] 2.824 0.066 2.694 2.824 2.959
beta0_pH[9,2] 3.950 0.212 3.554 3.940 4.394
beta0_pH[10,2] 3.994 0.125 3.755 3.990 4.252
beta0_pH[11,2] 1.885 0.009 1.868 1.885 1.903
beta0_pH[12,2] 1.793 0.067 1.663 1.793 1.926
beta0_pH[13,2] 2.496 0.026 2.445 2.496 2.549
beta0_pH[14,2] 1.804 0.024 1.759 1.804 1.850
beta0_pH[15,2] 2.545 0.017 2.511 2.545 2.579
beta0_pH[16,2] 2.766 0.025 2.716 2.765 2.814
beta0_pH[1,3] 1.106 0.365 0.208 1.169 1.663
beta0_pH[2,3] 2.127 0.017 2.094 2.126 2.161
beta0_pH[3,3] 2.692 0.024 2.645 2.692 2.740
beta0_pH[4,3] 2.980 0.035 2.914 2.980 3.048
beta0_pH[5,3] 3.325 2.205 0.836 2.786 8.983
beta0_pH[6,3] 1.349 0.370 0.697 1.324 2.165
beta0_pH[7,3] 0.861 0.050 0.764 0.860 0.957
beta0_pH[8,3] 0.248 0.294 -0.418 0.287 0.611
beta0_pH[9,3] -2.075 1.008 -4.362 -1.964 -0.423
beta0_pH[10,3] -0.411 0.794 -2.213 -0.318 0.898
beta0_pH[11,3] 1.112 0.526 0.073 1.148 2.035
beta0_pH[12,3] 0.949 0.654 -0.350 0.956 2.226
beta0_pH[13,3] 0.991 0.481 -0.009 1.021 1.897
beta0_pH[14,3] 0.562 0.446 -0.387 0.599 1.313
beta0_pH[15,3] 0.727 0.420 -0.126 0.748 1.463
beta0_pH[16,3] 1.174 0.629 -0.037 1.147 2.399
beta0_pH[11,4] 0.737 0.567 -0.765 0.800 1.633
beta0_pH[12,4] 0.722 0.668 -0.655 0.721 2.057
beta0_pH[13,4] 0.482 0.629 -1.128 0.585 1.450
beta0_pH[14,4] 0.302 0.469 -0.775 0.365 1.053
beta0_pH[15,4] 0.411 0.568 -1.088 0.517 1.239
beta0_pH[16,4] 0.708 0.535 -0.460 0.708 1.781
beta0_pH[11,5] 2.965 0.145 2.702 2.960 3.263
beta0_pH[12,5] -0.905 0.509 -1.945 -0.893 0.045
beta0_pH[13,5] 4.880 0.408 4.134 4.862 5.723
beta0_pH[14,5] 3.888 0.358 3.249 3.875 4.638
beta0_pH[15,5] 2.284 0.123 2.064 2.279 2.534
beta0_pH[16,5] 0.266 0.069 0.168 0.254 0.430
beta1_pH[1,1] 2.577 0.470 1.829 2.527 3.688
beta1_pH[2,1] 1.622 0.296 1.082 1.607 2.246
beta1_pH[3,1] 1.815 0.432 1.124 1.769 2.853
beta1_pH[4,1] 2.000 0.512 1.212 1.924 3.361
beta1_pH[5,1] 2.282 0.477 1.493 2.233 3.387
beta1_pH[6,1] 2.078 0.794 1.045 1.866 4.236
beta1_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,1] 2.324 0.877 1.176 2.090 4.614
beta1_pH[9,1] 1.972 0.578 1.190 1.860 3.619
beta1_pH[10,1] 2.113 0.341 1.508 2.085 2.853
beta1_pH[11,1] 3.315 0.541 2.368 3.258 4.531
beta1_pH[12,1] 2.928 0.285 2.376 2.927 3.493
beta1_pH[13,1] 3.452 0.249 3.001 3.441 3.979
beta1_pH[14,1] 2.845 0.322 2.236 2.840 3.502
beta1_pH[15,1] 2.410 0.288 1.898 2.393 3.022
beta1_pH[16,1] 4.006 0.303 3.419 3.999 4.622
beta1_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[1,3] 1.938 1.157 0.088 1.810 4.982
beta1_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,3] 4.760 2.489 0.964 4.210 11.108
beta1_pH[9,3] 7.176 2.574 3.281 6.756 13.251
beta1_pH[10,3] 6.308 2.515 2.686 5.852 12.012
beta1_pH[11,3] 2.181 0.575 1.073 2.167 3.313
beta1_pH[12,3] 3.266 0.677 2.020 3.240 4.661
beta1_pH[13,3] 2.000 0.516 0.991 1.978 3.070
beta1_pH[14,3] 2.229 0.501 1.348 2.197 3.271
beta1_pH[15,3] 2.097 0.484 1.212 2.075 3.077
beta1_pH[16,3] 3.702 0.682 2.355 3.728 5.003
beta1_pH[11,4] 2.071 0.613 1.008 2.030 3.536
beta1_pH[12,4] 3.321 0.684 1.979 3.304 4.770
beta1_pH[13,4] 1.526 0.638 0.486 1.434 3.090
beta1_pH[14,4] 2.086 0.536 1.156 2.021 3.266
beta1_pH[15,4] 1.610 0.599 0.667 1.531 3.122
beta1_pH[16,4] 3.757 0.623 2.476 3.768 5.106
beta1_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,1] 1.221 0.824 0.333 0.927 3.324
beta2_pH[2,1] 1.634 0.933 0.339 1.462 3.419
beta2_pH[3,1] 1.362 0.848 0.200 1.208 2.926
beta2_pH[4,1] 0.952 0.733 0.184 0.702 2.866
beta2_pH[5,1] 2.139 2.209 0.298 1.537 7.969
beta2_pH[6,1] 1.185 1.963 0.120 0.510 6.945
beta2_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,1] 0.976 1.646 0.144 0.414 6.007
beta2_pH[9,1] 1.179 1.676 0.173 0.630 6.278
beta2_pH[10,1] 1.434 1.762 0.320 0.855 6.979
beta2_pH[11,1] 0.459 0.337 0.179 0.369 1.214
beta2_pH[12,1] 1.553 0.820 0.619 1.285 3.658
beta2_pH[13,1] 1.022 0.428 0.548 0.915 2.199
beta2_pH[14,1] 1.117 0.670 0.425 0.917 3.182
beta2_pH[15,1] 1.045 0.626 0.380 0.876 2.933
beta2_pH[16,1] 1.062 0.338 0.631 1.003 1.875
beta2_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,3] 6.929 6.452 0.195 5.592 22.410
beta2_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,3] 1.530 1.374 0.142 1.086 5.505
beta2_pH[9,3] 0.303 0.567 0.091 0.208 0.985
beta2_pH[10,3] 0.455 0.738 0.115 0.276 2.119
beta2_pH[11,3] -5.744 4.524 -18.081 -4.431 -0.688
beta2_pH[12,3] -6.478 4.423 -18.063 -5.407 -1.141
beta2_pH[13,3] -5.823 4.392 -17.849 -4.546 -1.036
beta2_pH[14,3] -5.928 4.334 -18.249 -4.666 -1.227
beta2_pH[15,3] -6.135 4.433 -17.849 -4.904 -1.028
beta2_pH[16,3] -5.985 4.196 -17.328 -4.743 -1.446
beta2_pH[11,4] -6.064 4.578 -17.749 -4.877 -0.726
beta2_pH[12,4] -6.527 4.517 -17.763 -5.394 -1.125
beta2_pH[13,4] -6.197 4.548 -18.094 -5.022 -0.852
beta2_pH[14,4] -6.329 4.378 -17.532 -5.273 -1.108
beta2_pH[15,4] -6.327 4.515 -17.918 -5.181 -0.946
beta2_pH[16,4] -6.708 4.362 -18.349 -5.507 -1.562
beta2_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta3_pH[1,1] 34.914 1.087 32.844 34.926 37.206
beta3_pH[2,1] 32.502 1.376 29.900 32.489 35.235
beta3_pH[3,1] 32.898 1.518 29.593 32.907 36.146
beta3_pH[4,1] 32.763 1.968 29.565 32.582 37.426
beta3_pH[5,1] 27.659 1.167 26.018 27.493 30.989
beta3_pH[6,1] 32.971 3.380 25.497 33.045 40.749
beta3_pH[7,1] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,1] 35.871 3.830 29.192 36.199 42.551
beta3_pH[9,1] 31.965 2.624 27.786 31.618 39.100
beta3_pH[10,1] 32.487 1.233 30.155 32.454 35.077
beta3_pH[11,1] 31.370 1.109 29.318 31.382 33.498
beta3_pH[12,1] 30.775 0.465 29.822 30.801 31.654
beta3_pH[13,1] 31.412 0.535 30.429 31.398 32.526
beta3_pH[14,1] 32.564 0.758 31.285 32.514 34.147
beta3_pH[15,1] 30.943 0.753 29.416 30.957 32.389
beta3_pH[16,1] 30.130 0.371 29.395 30.138 30.853
beta3_pH[1,2] 34.430 5.516 25.539 34.394 43.549
beta3_pH[2,2] 34.571 5.444 25.508 34.712 43.486
beta3_pH[3,2] 30.511 7.200 18.570 30.461 42.338
beta3_pH[4,2] 30.473 7.182 18.586 30.489 42.269
beta3_pH[5,2] 30.422 7.209 18.638 30.307 42.415
beta3_pH[6,2] 30.489 7.150 18.607 30.487 42.354
beta3_pH[7,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[9,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[10,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[11,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[12,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[13,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[14,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[15,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[16,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[1,3] 38.383 3.265 30.669 39.730 42.717
beta3_pH[2,3] 34.609 5.453 25.570 34.713 43.599
beta3_pH[3,3] 30.397 7.173 18.672 30.180 42.267
beta3_pH[4,3] 30.443 7.234 18.540 30.433 42.319
beta3_pH[5,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[6,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[7,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,3] 42.028 1.342 38.391 42.059 43.908
beta3_pH[9,3] 36.154 1.920 34.529 35.427 41.722
beta3_pH[10,3] 37.842 2.475 34.650 37.167 43.347
beta3_pH[11,3] 40.709 1.048 38.941 40.698 42.657
beta3_pH[12,3] 42.036 0.702 40.366 42.264 42.835
beta3_pH[13,3] 41.402 0.709 40.008 41.420 42.782
beta3_pH[14,3] 40.460 0.665 39.138 40.459 41.734
beta3_pH[15,3] 40.633 0.658 39.122 40.652 41.925
beta3_pH[16,3] 41.856 0.388 41.036 41.907 42.621
beta3_pH[11,4] 40.704 1.034 38.980 40.598 42.996
beta3_pH[12,4] 41.971 0.740 40.344 42.210 42.904
beta3_pH[13,4] 42.027 1.072 39.858 42.116 43.805
beta3_pH[14,4] 40.491 0.681 39.075 40.508 41.754
beta3_pH[15,4] 41.375 1.012 39.527 41.277 43.390
beta3_pH[16,4] 41.536 0.368 40.923 41.522 42.287
beta3_pH[11,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[12,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[13,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[14,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[15,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[16,5] 25.000 0.000 25.000 25.000 25.000
beta4_pH[1,1] 0.263 0.021 0.223 0.263 0.304
beta4_pH[2,1] 0.211 0.005 0.201 0.211 0.221
beta4_pH[3,1] 0.173 0.007 0.158 0.173 0.186
beta4_pH[4,1] 0.294 0.011 0.271 0.294 0.316
beta4_pH[5,1] 0.420 0.269 0.084 0.389 0.912
beta4_pH[6,1] 0.775 0.468 0.333 0.724 1.447
beta4_pH[7,1] 0.868 0.348 0.362 0.809 1.692
beta4_pH[8,1] 0.159 0.069 0.045 0.150 0.309
beta4_pH[9,1] 0.214 0.086 0.062 0.209 0.394
beta4_pH[10,1] 0.001 0.001 0.000 0.001 0.003
beta4_pH[11,1] 0.447 0.012 0.425 0.447 0.471
beta4_pH[12,1] 0.641 0.022 0.598 0.640 0.685
beta4_pH[13,1] 0.631 0.013 0.606 0.632 0.658
beta4_pH[14,1] 0.670 0.030 0.613 0.669 0.731
beta4_pH[15,1] 0.797 0.021 0.755 0.797 0.838
beta4_pH[16,1] 1.221 0.042 1.144 1.219 1.309
beta4_pH[1,2] 0.306 0.094 0.140 0.300 0.508
beta4_pH[2,2] 0.513 0.022 0.470 0.512 0.557
beta4_pH[3,2] 0.446 0.012 0.423 0.445 0.469
beta4_pH[4,2] 0.654 0.038 0.580 0.654 0.732
beta4_pH[5,2] 0.733 0.883 0.109 0.540 2.561
beta4_pH[6,2] 0.933 0.928 0.197 0.745 2.667
beta4_pH[7,2] 0.888 0.402 0.354 0.816 1.849
beta4_pH[8,2] 0.844 0.203 0.469 0.837 1.271
beta4_pH[9,2] 1.164 1.131 0.308 0.873 3.842
beta4_pH[10,2] 0.456 0.050 0.358 0.457 0.557
beta4_pH[11,2] 1.270 0.051 1.173 1.269 1.372
beta4_pH[12,2] 0.580 0.108 0.386 0.573 0.805
beta4_pH[13,2] 0.943 0.043 0.864 0.943 1.029
beta4_pH[14,2] 0.929 0.063 0.808 0.929 1.056
beta4_pH[15,2] 1.011 0.025 0.961 1.011 1.061
beta4_pH[16,2] 1.405 0.075 1.263 1.405 1.556
beta4_pH[1,3] 0.286 0.088 0.130 0.282 0.462
beta4_pH[2,3] 0.052 0.017 0.021 0.051 0.086
beta4_pH[3,3] 0.052 0.010 0.033 0.052 0.071
beta4_pH[4,3] 0.172 0.032 0.110 0.172 0.236
beta4_pH[5,3] 0.744 0.708 0.123 0.562 2.654
beta4_pH[6,3] 0.902 0.967 0.203 0.717 2.710
beta4_pH[7,3] 0.900 0.456 0.358 0.818 1.981
beta4_pH[8,3] 0.420 0.199 0.113 0.394 0.886
beta4_pH[9,3] 0.943 0.942 0.173 0.707 3.035
beta4_pH[10,3] 0.061 0.041 0.003 0.055 0.154
beta4_pH[11,3] 0.363 0.016 0.333 0.363 0.394
beta4_pH[12,3] 0.440 0.042 0.360 0.439 0.527
beta4_pH[13,3] 0.838 0.037 0.765 0.838 0.914
beta4_pH[14,3] 0.823 0.066 0.700 0.821 0.960
beta4_pH[15,3] 0.635 0.017 0.603 0.635 0.668
beta4_pH[16,3] 0.637 0.036 0.569 0.635 0.708
beta4_pH[1,4] 0.628 0.328 0.209 0.563 1.412
beta4_pH[2,4] 0.802 0.792 0.113 0.588 2.835
beta4_pH[3,4] 0.846 0.880 0.113 0.594 3.148
beta4_pH[4,4] 0.732 0.474 0.218 0.629 1.814
beta4_pH[5,4] 0.770 0.826 0.118 0.559 2.763
beta4_pH[6,4] 0.946 1.028 0.215 0.736 2.914
beta4_pH[7,4] 0.881 0.401 0.353 0.815 1.821
beta4_pH[8,4] 0.723 0.587 0.159 0.594 2.096
beta4_pH[9,4] 0.795 0.949 0.122 0.576 2.743
beta4_pH[10,4] 1.519 3.112 0.006 0.392 11.974
beta4_pH[11,4] 0.444 0.017 0.411 0.444 0.478
beta4_pH[12,4] 0.555 0.053 0.457 0.554 0.663
beta4_pH[13,4] 0.925 0.036 0.853 0.925 0.997
beta4_pH[14,4] 0.917 0.075 0.778 0.915 1.074
beta4_pH[15,4] 0.714 0.020 0.676 0.713 0.755
beta4_pH[16,4] 0.685 0.041 0.609 0.683 0.766
beta4_pH[1,5] 0.629 0.332 0.206 0.564 1.423
beta4_pH[2,5] 0.826 1.026 0.114 0.593 2.852
beta4_pH[3,5] 0.875 0.984 0.115 0.612 3.501
beta4_pH[4,5] 0.734 0.465 0.209 0.628 1.914
beta4_pH[5,5] 0.737 0.759 0.104 0.566 2.372
beta4_pH[6,5] 0.908 0.716 0.216 0.737 2.708
beta4_pH[7,5] 0.894 0.505 0.346 0.813 1.893
beta4_pH[8,5] 0.751 0.743 0.151 0.590 2.299
beta4_pH[9,5] 0.778 0.881 0.122 0.579 2.731
beta4_pH[10,5] 1.536 3.217 0.007 0.392 12.349
beta4_pH[11,5] 0.603 0.051 0.508 0.601 0.709
beta4_pH[12,5] 0.587 0.171 0.307 0.566 0.968
beta4_pH[13,5] 0.776 0.062 0.661 0.774 0.903
beta4_pH[14,5] 0.840 0.115 0.631 0.832 1.078
beta4_pH[15,5] 1.144 0.073 1.011 1.144 1.296
beta4_pH[16,5] 17.919 4.039 9.965 18.068 24.603
beta0_pelagic[1] 2.266 0.036 2.199 2.266 2.341
beta0_pelagic[2] 1.538 0.007 1.524 1.538 1.552
beta0_pelagic[3] 0.210 0.262 -0.536 0.272 0.569
beta0_pelagic[4] 0.232 0.217 -0.292 0.249 0.605
beta0_pelagic[5] 1.452 0.156 1.154 1.451 1.762
beta0_pelagic[6] 1.353 0.059 1.240 1.352 1.472
beta0_pelagic[7] 1.483 0.038 1.413 1.482 1.559
beta0_pelagic[8] 1.673 0.042 1.594 1.673 1.759
beta0_pelagic[9] 2.795 0.065 2.671 2.794 2.926
beta0_pelagic[10] 2.603 0.031 2.545 2.603 2.667
beta0_pelagic[11] -0.165 0.160 -0.542 -0.148 0.105
beta0_pelagic[12] 1.802 0.035 1.735 1.801 1.874
beta0_pelagic[13] -0.833 0.464 -1.682 -0.812 -0.001
beta0_pelagic[14] -0.275 0.170 -0.641 -0.263 0.027
beta0_pelagic[15] -0.208 0.034 -0.299 -0.195 -0.181
beta0_pelagic[16] -0.584 0.393 -1.463 -0.525 -0.007
beta1_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[3] 0.735 0.284 0.305 0.692 1.544
beta1_pelagic[4] 0.815 0.246 0.391 0.798 1.414
beta1_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[11] 1.181 0.224 0.810 1.161 1.705
beta1_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[13] 2.596 0.529 1.454 2.629 3.445
beta1_pelagic[14] 1.758 0.416 1.117 1.684 2.800
beta1_pelagic[15] 0.079 0.146 0.000 0.012 0.566
beta1_pelagic[16] 1.777 0.581 0.867 1.713 2.953
beta2_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[3] 6.596 6.088 0.223 5.027 21.132
beta2_pelagic[4] 7.696 5.929 0.366 6.492 21.586
beta2_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[11] 1.799 1.585 0.294 1.249 6.237
beta2_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[13] 0.196 0.090 0.106 0.173 0.424
beta2_pelagic[14] 0.882 0.854 0.221 0.605 3.387
beta2_pelagic[15] 2.643 1.882 0.174 2.276 7.063
beta2_pelagic[16] 0.565 0.837 0.115 0.286 3.366
beta3_pelagic[1] 28.747 3.842 22.396 28.655 35.109
beta3_pelagic[2] 28.741 3.906 22.294 28.786 35.128
beta3_pelagic[3] 27.940 2.695 23.138 28.383 32.104
beta3_pelagic[4] 26.099 1.214 23.217 26.304 28.468
beta3_pelagic[5] 29.459 6.750 18.673 29.217 40.536
beta3_pelagic[6] 29.696 6.657 18.584 29.859 40.515
beta3_pelagic[7] 29.652 6.729 18.558 29.693 40.519
beta3_pelagic[8] 29.395 6.589 18.477 29.383 40.302
beta3_pelagic[9] 29.631 6.627 18.534 29.718 40.449
beta3_pelagic[10] 29.345 6.675 18.500 29.244 40.449
beta3_pelagic[11] 30.015 1.185 27.793 29.947 32.491
beta3_pelagic[12] 30.000 0.000 30.000 30.000 30.000
beta3_pelagic[13] 31.843 4.165 24.861 31.539 40.019
beta3_pelagic[14] 34.174 1.587 31.558 33.987 37.877
beta3_pelagic[15] 33.959 5.636 23.728 34.685 41.684
beta3_pelagic[16] 31.009 2.907 25.298 30.946 37.453
beta0_yellow[1] -0.482 0.080 -0.636 -0.483 -0.328
beta0_yellow[2] 0.540 0.084 0.348 0.548 0.683
beta0_yellow[3] -0.208 0.129 -0.457 -0.213 0.016
beta0_yellow[4] 0.889 0.264 0.104 0.964 1.190
beta0_yellow[5] 0.254 0.324 -0.439 0.282 0.834
beta0_yellow[6] 0.819 0.106 0.618 0.816 1.029
beta0_yellow[7] 0.998 0.094 0.818 0.997 1.184
beta0_yellow[8] 1.026 0.094 0.845 1.024 1.212
beta0_yellow[9] 0.461 0.098 0.266 0.463 0.652
beta0_yellow[10] 0.536 0.068 0.402 0.537 0.673
beta0_yellow[11] -0.198 0.128 -0.485 -0.188 0.027
beta0_yellow[12] -1.081 0.092 -1.266 -1.079 -0.901
beta0_yellow[13] -1.352 0.169 -1.732 -1.344 -1.059
beta0_yellow[14] -0.251 0.153 -0.585 -0.232 -0.003
beta0_yellow[15] -1.250 0.134 -1.569 -1.234 -1.031
beta0_yellow[16] -0.252 0.093 -0.437 -0.252 -0.071
beta1_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[2] 1.010 0.235 0.562 1.006 1.477
beta1_yellow[3] 0.465 0.246 0.019 0.478 0.929
beta1_yellow[4] 1.321 0.598 0.624 1.170 3.206
beta1_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[11] 0.309 0.395 0.001 0.117 1.267
beta1_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[13] 1.101 1.115 0.022 0.785 4.577
beta1_yellow[14] 0.718 1.018 0.013 0.423 4.151
beta1_yellow[15] 0.163 0.200 0.011 0.080 0.680
beta1_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[2] -7.472 5.703 -21.226 -6.271 -0.470
beta2_yellow[3] -7.870 5.872 -21.917 -6.543 -0.444
beta2_yellow[4] -5.748 5.901 -20.255 -3.806 -0.119
beta2_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[11] -7.804 5.210 -20.927 -7.106 -0.608
beta2_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[13] -5.996 5.391 -19.862 -4.438 -0.311
beta2_yellow[14] -7.048 5.580 -20.745 -5.857 -0.382
beta2_yellow[15] -6.996 5.049 -18.958 -6.075 -0.427
beta2_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta3_yellow[1] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[2] 28.938 1.190 27.376 28.700 32.365
beta3_yellow[3] 32.373 3.372 23.182 32.586 38.686
beta3_yellow[4] 30.026 3.324 24.088 28.844 35.889
beta3_yellow[5] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[6] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[7] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[8] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[9] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[10] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[11] 34.514 3.421 29.688 33.679 42.153
beta3_yellow[12] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[13] 33.622 2.755 29.263 33.671 40.182
beta3_yellow[14] 33.978 3.649 29.189 33.452 43.035
beta3_yellow[15] 36.492 3.833 30.374 36.341 43.469
beta3_yellow[16] 25.000 0.000 25.000 25.000 25.000
beta0_black[1] 0.009 0.075 -0.137 0.007 0.151
beta0_black[2] 1.854 0.032 1.791 1.854 1.918
beta0_black[3] 1.240 0.037 1.168 1.240 1.312
beta0_black[4] 2.423 0.047 2.334 2.423 2.519
beta0_black[5] 5.262 2.487 1.850 4.797 11.496
beta0_black[6] 5.292 2.169 2.670 4.748 10.918
beta0_black[7] 4.094 2.073 1.578 3.557 9.709
beta0_black[8] 0.882 0.216 0.471 0.880 1.331
beta0_black[9] 2.816 0.276 2.300 2.810 3.377
beta0_black[10] 1.285 0.051 1.187 1.285 1.387
beta0_black[11] 2.977 0.051 2.877 2.977 3.078
beta0_black[12] 4.917 0.187 4.553 4.919 5.278
beta0_black[13] 0.438 0.335 0.017 0.357 1.225
beta0_black[14] 2.853 0.104 2.656 2.853 3.056
beta0_black[15] 0.722 0.052 0.620 0.722 0.822
beta0_black[16] 3.773 0.073 3.628 3.771 3.920
beta2_black[1] 0.000 0.000 0.000 0.000 0.000
beta2_black[2] 0.000 0.000 0.000 0.000 0.000
beta2_black[3] 0.000 0.000 0.000 0.000 0.000
beta2_black[4] 0.000 0.000 0.000 0.000 0.000
beta2_black[5] 0.000 0.000 0.000 0.000 0.000
beta2_black[6] 0.000 0.000 0.000 0.000 0.000
beta2_black[7] 0.000 0.000 0.000 0.000 0.000
beta2_black[8] 0.000 0.000 0.000 0.000 0.000
beta2_black[9] 0.000 0.000 0.000 0.000 0.000
beta2_black[10] 0.000 0.000 0.000 0.000 0.000
beta2_black[11] 0.000 0.000 0.000 0.000 0.000
beta2_black[12] 0.000 0.000 0.000 0.000 0.000
beta2_black[13] -2.514 1.798 -7.071 -2.078 -0.285
beta2_black[14] 0.000 0.000 0.000 0.000 0.000
beta2_black[15] 0.000 0.000 0.000 0.000 0.000
beta2_black[16] 0.000 0.000 0.000 0.000 0.000
beta3_black[1] 25.000 0.000 25.000 25.000 25.000
beta3_black[2] 25.000 0.000 25.000 25.000 25.000
beta3_black[3] 25.000 0.000 25.000 25.000 25.000
beta3_black[4] 25.000 0.000 25.000 25.000 25.000
beta3_black[5] 25.000 0.000 25.000 25.000 25.000
beta3_black[6] 25.000 0.000 25.000 25.000 25.000
beta3_black[7] 25.000 0.000 25.000 25.000 25.000
beta3_black[8] 25.000 0.000 25.000 25.000 25.000
beta3_black[9] 25.000 0.000 25.000 25.000 25.000
beta3_black[10] 25.000 0.000 25.000 25.000 25.000
beta3_black[11] 25.000 0.000 25.000 25.000 25.000
beta3_black[12] 25.000 0.000 25.000 25.000 25.000
beta3_black[13] 37.345 3.307 27.209 38.125 40.990
beta3_black[14] 25.000 0.000 25.000 25.000 25.000
beta3_black[15] 25.000 0.000 25.000 25.000 25.000
beta3_black[16] 25.000 0.000 25.000 25.000 25.000
beta4_black[1] -0.325 0.091 -0.501 -0.324 -0.149
beta4_black[2] 0.257 0.052 0.153 0.257 0.361
beta4_black[3] -0.719 0.063 -0.844 -0.718 -0.594
beta4_black[4] 0.626 0.132 0.373 0.623 0.892
beta4_black[5] 0.753 1.450 -1.308 0.479 4.491
beta4_black[6] 0.857 1.543 -0.981 0.528 4.584
beta4_black[7] 0.815 1.519 -0.969 0.518 4.227
beta4_black[8] 0.087 0.282 -0.509 0.095 0.602
beta4_black[9] 1.048 0.838 -0.107 0.875 3.206
beta4_black[10] 0.075 0.066 -0.054 0.075 0.203
beta4_black[11] -0.608 0.065 -0.732 -0.610 -0.478
beta4_black[12] 1.230 0.288 0.679 1.222 1.824
beta4_black[13] -1.002 0.077 -1.158 -0.999 -0.853
beta4_black[14] 0.219 0.178 -0.119 0.215 0.578
beta4_black[15] 0.002 0.046 -0.089 0.003 0.093
beta4_black[16] -0.412 0.117 -0.644 -0.411 -0.183
beta0_dsr[11] -0.329 0.010 -0.347 -0.328 -0.310
beta0_dsr[12] 1.409 0.044 1.322 1.410 1.492
beta0_dsr[13] 0.155 0.021 0.113 0.155 0.196
beta0_dsr[14] -0.127 0.017 -0.161 -0.127 -0.094
beta0_dsr[15] 0.205 0.013 0.179 0.205 0.231
beta0_dsr[16] 0.165 0.016 0.133 0.164 0.196
beta1_dsr[11] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[12] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[13] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[14] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[15] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[16] 0.000 0.000 0.000 0.000 0.000
beta4_dsr[11] 0.544 0.035 0.476 0.544 0.613
beta4_dsr[12] 0.794 0.223 0.382 0.789 1.254
beta4_dsr[13] -0.078 0.040 -0.156 -0.078 0.000
beta4_dsr[14] -0.206 0.071 -0.349 -0.207 -0.066
beta4_dsr[15] 0.216 0.021 0.175 0.216 0.258
beta4_dsr[16] -0.055 0.045 -0.143 -0.054 0.035
beta0_slope[11] -2.209 0.021 -2.252 -2.209 -2.168
beta0_slope[12] -4.987 0.187 -5.377 -4.985 -4.635
beta0_slope[13] -1.502 0.034 -1.568 -1.501 -1.434
beta0_slope[14] -2.712 0.050 -2.808 -2.711 -2.615
beta0_slope[15] -1.592 0.022 -1.635 -1.593 -1.550
beta0_slope[16] -2.657 0.038 -2.731 -2.657 -2.581
beta1_slope[11] 0.000 0.000 0.000 0.000 0.000
beta1_slope[12] 0.000 0.000 0.000 0.000 0.000
beta1_slope[13] 0.000 0.000 0.000 0.000 0.000
beta1_slope[14] 0.000 0.000 0.000 0.000 0.000
beta1_slope[15] 0.000 0.000 0.000 0.000 0.000
beta1_slope[16] 0.000 0.000 0.000 0.000 0.000
beta4_slope[11] -0.333 0.074 -0.484 -0.332 -0.189
beta4_slope[12] -4.684 0.789 -6.345 -4.637 -3.252
beta4_slope[13] 0.328 0.053 0.225 0.329 0.433
beta4_slope[14] -0.159 0.151 -0.466 -0.155 0.129
beta4_slope[15] -0.359 0.033 -0.424 -0.359 -0.295
beta4_slope[16] -0.137 0.102 -0.329 -0.136 0.063
beta5_slope[11] 3.065 0.108 2.860 3.063 3.289
beta5_slope[12] 2.006 0.388 1.279 1.994 2.807
beta5_slope[13] 0.879 0.080 0.721 0.877 1.034
beta5_slope[14] 3.743 0.354 3.091 3.730 4.471
beta5_slope[15] 1.308 0.073 1.173 1.307 1.456
beta5_slope[16] 2.718 0.236 2.288 2.713 3.192
tau_comp[1] 3.808 0.294 3.264 3.801 4.407
tau_comp[2] 2.857 0.253 2.376 2.846 3.381
tau_comp[3] 0.873 0.074 0.740 0.870 1.026
tau_comp[4] 1.782 0.273 1.315 1.761 2.348
tau_comp[5] 1.288 0.187 0.956 1.277 1.683
sigma_H[1] 0.207 0.049 0.123 0.205 0.314
sigma_H[2] 0.174 0.029 0.127 0.171 0.239
sigma_H[3] 0.187 0.040 0.115 0.184 0.274
sigma_H[4] 0.461 0.079 0.325 0.453 0.643
sigma_H[5] 0.958 0.201 0.591 0.949 1.375
sigma_H[6] 0.537 0.199 0.139 0.531 0.948
sigma_H[7] 0.311 0.063 0.212 0.302 0.453
sigma_H[8] 0.421 0.089 0.284 0.408 0.613
sigma_H[9] 0.510 0.114 0.330 0.494 0.762
sigma_H[10] 0.212 0.038 0.147 0.208 0.297
sigma_H[11] 0.265 0.043 0.191 0.261 0.361
sigma_H[12] 0.380 0.147 0.200 0.335 0.729
sigma_H[13] 0.218 0.036 0.157 0.215 0.297
sigma_H[14] 0.497 0.087 0.347 0.491 0.682
sigma_H[15] 0.230 0.037 0.166 0.227 0.310
sigma_H[16] 0.215 0.039 0.148 0.210 0.304
lambda_H[1] 2.187 3.432 0.075 1.003 11.711
lambda_H[2] 5.324 5.074 0.498 3.827 19.283
lambda_H[3] 5.476 8.207 0.300 2.848 27.386
lambda_H[4] 0.006 0.005 0.001 0.005 0.019
lambda_H[5] 2.805 7.399 0.029 0.660 21.269
lambda_H[6] 6.698 14.266 0.008 0.439 48.453
lambda_H[7] 0.012 0.009 0.002 0.010 0.033
lambda_H[8] 7.391 9.914 0.041 3.842 34.176
lambda_H[9] 0.015 0.010 0.003 0.013 0.041
lambda_H[10] 0.233 0.447 0.023 0.154 0.848
lambda_H[11] 0.223 0.350 0.010 0.099 1.111
lambda_H[12] 3.766 5.010 0.169 2.133 17.755
lambda_H[13] 3.165 2.932 0.263 2.315 10.782
lambda_H[14] 3.312 3.933 0.204 2.172 13.004
lambda_H[15] 0.018 0.019 0.003 0.014 0.061
lambda_H[16] 0.551 0.953 0.033 0.265 2.923
mu_lambda_H[1] 4.126 1.925 1.066 3.925 8.362
mu_lambda_H[2] 3.627 1.951 0.533 3.437 7.783
mu_lambda_H[3] 3.350 1.845 0.726 3.030 7.480
sigma_lambda_H[1] 8.260 4.435 1.734 7.602 18.243
sigma_lambda_H[2] 8.012 4.711 0.819 7.412 18.064
sigma_lambda_H[3] 6.155 4.031 0.958 5.144 16.043
beta_H[1,1] 6.475 1.371 3.089 6.702 8.465
beta_H[2,1] 9.753 0.611 8.459 9.788 10.818
beta_H[3,1] 7.853 0.753 6.111 7.924 9.103
beta_H[4,1] 10.665 7.648 -5.155 10.832 25.259
beta_H[5,1] -0.269 2.600 -5.605 -0.086 4.219
beta_H[6,1] 3.308 4.240 -7.235 4.864 8.354
beta_H[7,1] 1.170 5.691 -11.193 1.487 11.506
beta_H[8,1] 1.440 4.332 -2.441 1.182 3.865
beta_H[9,1] 13.660 5.531 2.909 13.496 24.939
beta_H[10,1] 6.796 1.814 3.195 6.877 10.353
beta_H[11,1] 4.828 3.541 -2.877 5.351 10.003
beta_H[12,1] 2.583 1.087 0.728 2.511 4.855
beta_H[13,1] 9.098 0.879 7.230 9.170 10.542
beta_H[14,1] 2.131 1.004 0.113 2.121 4.147
beta_H[15,1] -7.046 3.319 -13.180 -7.181 0.009
beta_H[16,1] 3.849 2.899 -0.771 3.510 10.689
beta_H[1,2] 7.778 0.276 7.231 7.791 8.298
beta_H[2,2] 9.888 0.169 9.545 9.889 10.207
beta_H[3,2] 8.802 0.198 8.400 8.800 9.190
beta_H[4,2] 3.035 1.407 0.367 2.999 5.883
beta_H[5,2] 1.745 0.989 -0.239 1.713 3.615
beta_H[6,2] 6.042 1.101 3.538 6.224 7.765
beta_H[7,2] 2.427 1.070 0.488 2.381 4.726
beta_H[8,2] 2.945 1.207 0.887 3.121 4.220
beta_H[9,2] 3.134 1.073 1.089 3.128 5.262
beta_H[10,2] 8.028 0.373 7.215 8.045 8.729
beta_H[11,2] 9.834 0.640 8.835 9.740 11.224
beta_H[12,2] 3.929 0.384 3.202 3.918 4.693
beta_H[13,2] 9.153 0.241 8.713 9.148 9.630
beta_H[14,2] 4.016 0.336 3.371 4.017 4.704
beta_H[15,2] 11.546 0.597 10.318 11.568 12.652
beta_H[16,2] 4.477 0.802 3.011 4.459 6.130
beta_H[1,3] 8.457 0.267 7.970 8.445 8.999
beta_H[2,3] 9.873 0.145 9.588 9.871 10.173
beta_H[3,3] 9.470 0.164 9.160 9.470 9.805
beta_H[4,3] -2.170 0.863 -3.940 -2.149 -0.530
beta_H[5,3] 3.878 0.591 2.677 3.880 5.010
beta_H[6,3] 8.334 1.202 6.687 7.969 10.949
beta_H[7,3] -2.587 0.642 -3.943 -2.581 -1.357
beta_H[8,3] 5.335 0.590 4.687 5.235 6.684
beta_H[9,3] -2.665 0.689 -4.040 -2.653 -1.392
beta_H[10,3] 8.563 0.276 8.046 8.553 9.133
beta_H[11,3] 8.503 0.295 7.876 8.523 9.018
beta_H[12,3] 5.269 0.327 4.541 5.312 5.828
beta_H[13,3] 8.861 0.175 8.514 8.870 9.198
beta_H[14,3] 5.785 0.267 5.195 5.806 6.257
beta_H[15,3] 10.282 0.288 9.742 10.279 10.871
beta_H[16,3] 5.965 0.610 4.698 5.994 7.077
beta_H[1,4] 8.097 0.186 7.713 8.100 8.443
beta_H[2,4] 9.861 0.128 9.593 9.869 10.103
beta_H[3,4] 9.971 0.154 9.640 9.979 10.244
beta_H[4,4] 11.812 0.471 10.872 11.810 12.722
beta_H[5,4] 5.581 0.734 4.398 5.501 7.249
beta_H[6,4] 7.205 0.904 5.309 7.430 8.488
beta_H[7,4] 8.549 0.354 7.839 8.558 9.231
beta_H[8,4] 6.707 0.275 6.135 6.728 7.129
beta_H[9,4] 7.420 0.448 6.561 7.410 8.356
beta_H[10,4] 7.524 0.230 7.075 7.514 7.988
beta_H[11,4] 9.504 0.195 9.139 9.502 9.915
beta_H[12,4] 7.249 0.198 6.866 7.247 7.664
beta_H[13,4] 9.070 0.138 8.796 9.074 9.328
beta_H[14,4] 7.837 0.212 7.426 7.832 8.256
beta_H[15,4] 9.468 0.217 9.059 9.469 9.895
beta_H[16,4] 9.419 0.243 8.967 9.414 9.892
beta_H[1,5] 8.896 0.158 8.567 8.902 9.184
beta_H[2,5] 10.678 0.100 10.490 10.673 10.886
beta_H[3,5] 10.767 0.159 10.479 10.762 11.089
beta_H[4,5] 8.351 0.501 7.386 8.343 9.385
beta_H[5,5] 5.193 0.608 3.695 5.274 6.193
beta_H[6,5] 8.928 0.637 7.971 8.801 10.371
beta_H[7,5] 6.571 0.351 5.893 6.563 7.300
beta_H[8,5] 8.250 0.229 7.907 8.226 8.710
beta_H[9,5] 8.195 0.442 7.307 8.193 9.060
beta_H[10,5] 10.035 0.227 9.591 10.032 10.474
beta_H[11,5] 11.602 0.226 11.166 11.598 12.067
beta_H[12,5] 8.604 0.186 8.222 8.602 8.971
beta_H[13,5] 10.141 0.133 9.888 10.136 10.409
beta_H[14,5] 9.242 0.218 8.836 9.236 9.688
beta_H[15,5] 11.317 0.232 10.854 11.322 11.763
beta_H[16,5] 9.838 0.188 9.466 9.842 10.192
beta_H[1,6] 10.127 0.247 9.738 10.087 10.726
beta_H[2,6] 11.469 0.113 11.246 11.468 11.693
beta_H[3,6] 10.615 0.159 10.273 10.623 10.912
beta_H[4,6] 12.928 0.881 11.176 12.951 14.572
beta_H[5,6] 5.673 0.627 4.539 5.648 6.990
beta_H[6,6] 8.721 0.559 7.446 8.755 9.736
beta_H[7,6] 10.560 0.599 9.362 10.577 11.710
beta_H[8,6] 9.517 0.287 8.976 9.538 9.978
beta_H[9,6] 8.268 0.751 6.823 8.257 9.830
beta_H[10,6] 9.150 0.342 8.419 9.174 9.740
beta_H[11,6] 10.687 0.368 9.903 10.718 11.333
beta_H[12,6] 9.375 0.237 8.954 9.358 9.870
beta_H[13,6] 11.094 0.158 10.794 11.087 11.425
beta_H[14,6] 9.916 0.278 9.333 9.924 10.433
beta_H[15,6] 10.551 0.414 9.743 10.552 11.371
beta_H[16,6] 10.598 0.259 10.045 10.610 11.094
beta_H[1,7] 10.113 1.282 6.903 10.423 11.824
beta_H[2,7] 12.108 0.513 10.979 12.135 13.061
beta_H[3,7] 10.268 0.672 8.782 10.339 11.371
beta_H[4,7] 1.305 4.499 -7.362 1.279 10.626
beta_H[5,7] 5.871 1.988 1.997 5.933 9.845
beta_H[6,7] 7.198 3.182 -0.944 8.362 11.039
beta_H[7,7] 7.634 3.046 1.705 7.650 13.745
beta_H[8,7] 10.911 1.102 9.155 10.894 12.749
beta_H[9,7] 5.070 3.900 -2.614 5.051 12.742
beta_H[10,7] 9.992 1.621 7.062 9.847 13.489
beta_H[11,7] 11.376 1.807 8.237 11.174 15.393
beta_H[12,7] 9.776 0.986 7.553 9.857 11.449
beta_H[13,7] 11.721 0.720 10.048 11.794 12.920
beta_H[14,7] 10.646 0.915 8.732 10.668 12.368
beta_H[15,7] 13.949 2.183 9.687 13.957 18.318
beta_H[16,7] 12.807 1.394 10.359 12.654 15.906
beta0_H[1] 9.004 18.415 -27.598 8.992 46.796
beta0_H[2] 10.274 7.997 -6.558 10.199 26.334
beta0_H[3] 9.248 9.611 -11.443 9.320 28.395
beta0_H[4] 6.726 183.633 -364.760 10.190 361.570
beta0_H[5] 3.367 27.012 -53.826 4.186 56.229
beta0_H[6] 7.466 51.416 -110.292 7.744 122.317
beta0_H[7] 4.313 141.035 -261.517 4.743 290.653
beta0_H[8] 6.129 31.125 -22.341 6.327 35.134
beta0_H[9] 8.763 117.271 -223.323 7.429 244.529
beta0_H[10] 8.999 37.464 -65.352 8.742 86.207
beta0_H[11] 9.892 52.166 -93.342 10.149 121.420
beta0_H[12] 6.415 12.797 -18.955 6.672 29.703
beta0_H[13] 9.791 13.048 -11.921 9.875 32.673
beta0_H[14] 7.120 11.743 -16.856 7.012 31.803
beta0_H[15] 9.560 113.440 -219.703 11.090 228.993
beta0_H[16] 8.272 30.581 -53.495 7.719 70.437